Calculation of the interaction of a neutron spin with an atomic electric field 

R. Golub^, S.K. Lamoreaux^ 
1. Hahn-Meitner Institute, Glienicker Str. 100, D 1^109 Berlin, 

Germany 

2. Los Alamos National Laboratory, Physics Division, P-23, 

M.S. H803, Los Alamos, NM 87545 

(February 8, 2008) 



0^ 



o 

in 



The Thomas- Fermi approximation for an atomic wavefunc- 
tion is used to calculate the interaction of a neutron spin with 
the atomic electric field, either through the motional magnetic 
(« X E) or possibly electric (due to the possible existence of 
a neutron permanent electric dipole moment) couplings. 



I. INTRODUCTION 



Because the calculations regarding this effect do not, 
to our knowledge, appear in the literature and the cal- 
culations regarding the SchuU experiment are lacking in 
detail, we felt it worthwhile to estimate the size of a neu- 
tron spin rotation due to an EDM interaction, and in 
addition, include the analysis of the v x E interaction 
originally studied by Schwinger ^ and demonstrated by 
ShuU and Nathans. [g| We use the Thomas-Fermi model 
of the atom to give the approximate atomic electric field. 



II. THOMAS-FERMI MODEL OF THE ATOM 
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Most recent experiments to search for the neutron elec- 
tric dipole moment (EDM) involve the neutron interact- 
ing with electric fields created by laboratory apparatus. 
However there is also the possibility of using electric fields 
produced by atoms in crystals. There are some reasons 
for believing this might be advantageous - atomic fields 
are large and the neutrons interact with many atoms in 
a coherent fashion. However as we will show the mea- 
surable effects are quite small compared to known back- 
ground effects due to the motional electric field. 

The effect in a crystal experiment can be enhanced 
if the scatteing amplitude due to the EDM interaction 
can be made to interfere with the much larger scattering 
amplitude of a nucleus. It will be seen below that the 
scattering amplitude due to the EDM interaction is imag- 
inary so it can only interfere with an imaginary nuclear 
amplitude. The idea of searching for a neutron EDM by 
measuring the interference between the scattering from 
an atomic electric field (due to the EDM interaction) with 
nuclear scattering was proposed by ShuU; an experiment 
was performed in 1964 and was based on scattering from 
a CdS crystal, because Cd, a strong absorber, has a large 
imaginary scattering amplitude. M In this experiment, 
the penetration depth, hence reflectivity, depends on the 
orientation of the neutron spin relative to the momentum 
transfer. ||J| 

Recently, a proposal to search for a neutron EDM by 
scattering in a perfect Si crystal was put forward. [^ In 
this case, the imaginary part of the nuclear scattering 
length is very small, and the proposed observable is a 
rotation of the neutron spin direction caused by the su- 
perpostion of the spin dependent imaginary amplitude 
with the real nuclear amplitude. 



The Thomas-Fermi model of the atom is fully devel- 
oped in §70 of M. Briefly, the atom is treated semi- 
classically with the electron density as a function of po- 
sition determined by phase space considerations. This 
leads to a universal function (i.e., it does not depend 
on atomic number Z) for the self-consistent electric field 
within the atom, 



d\/dx'^x'^' 



(1) 



where x describes the shielding (assumed spherically 
symmetric) of the nuclear point charge, with the bound- 
ary conditions that x(0) — 1 and x(oo) = (the latter 
condition determins x'(0)), and the radius r is related to 
X as 



= xbZ 



-1/3. 
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(where we are using atomic units so m^e^ /h 
electric potential within an atom is given by 

,, , Ze frZ^/^\ Z'^'^xix) 
cj){r) = — X 
r 
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We point out that the Thomas-Fermi model does not 
apply for either very large or very small x; however, the 
major contribution to the scattering integral is from x « 
1, and we would expect this model to give reasonably 
accurate results. 



III. INTERACTION OF AN EDM WITH THE 
ATOMIC ELECTRIC FIELD 

We are interested in calculating the spin-dependent 
neutron scattering length by the Born approximation for 



either the v x E field or a neutron permanent electric 
dipole moment (EDM). Consider first the EDM interac- 
tion, 



V{r) = -dea ■ E{r) 



(4) 



where d is the dipole moment length, e is the magnitude 
of the electron charge, a is a Pauli matrix, and E{r) is 
an electric field. The scattering amplitude (length) can 
be determined by use of the Born approximation. 
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IV. y X £ INTERACTION 

Next, consider the vxE motional magnetic field inter- 
action which couples to the neutron magnetic moment, 
which was first considered by Schwinger in 1948. H The 
possibility of measuring effects from the motional field 
has been discussed in regard to non-centrosymmetric cry- 
tals (a— quartz) in which case a non-zero average electric 
field between scattering planes can exist. H,g] However, 
as has been pointed out, there is avx E observable even 
for symmetric crystals. [^ 

The hamiltonian for the v x E interaction is 



Taking the momentum transfer q along z as the quanti- 
zation axis, and using the fact that the electric field is 
spherically symmetric, we find 

/•oo 

a^cr,^ deE{r)e'''''''°''^ cos sine r^drde (6) 
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and the other components are zero because of symmetry. 
Taking E(r) = —(f>'(r)f, and using the Thomas-Fermi 
wavefunction, the scattering length can be written as, 
taking f3 == bZ-^/^, 
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For the Born approximation, we use the matrix element 
{k2\V{f)\ki) = M . i(fc2|^ X E{f) ~ E{r) X ^\h) 

2 ninC TTlnC 

(12) 



where ki, fc2 label the incoming and outgoing neutron 
wavefunction. Taking 
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which can be rewritten as 
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where f{j3q) is the imaginary part (the real part is zero) 
of the dimcnsionlcss integral in the previous equation. 
The Thomas-Fermi equation was numerically solved us- 
ing a Runge-Kutta technique, and the integral numer- 
ically evaluated. The results, as a function of fiq, are 
shown in Fig. 1. 

For the case of Si, Z — lA, (3 — 0.367; a typical q is 
approximately 27r/2A x 0.5A/a.u. giving Pq = 0.58, and 
the dimensionless integral is about 1. Thus, the difference 
in the scattering length for the two spin states (along ztq) 
is 



Again, assume that g lies along z; we note that q 
is perpendicular to fci + fc2 because 



(fc2 - fcl) • (/Cl + k2) 



^2 



kl 







(13) 
(14) 

(15) 
k2-ki 

(16) 



for elastic scattering. By symmetry, the effective electric 
field lies along z, and if we assume ki + k2 is along y 
and q (and has magnitude 2k cos Og , 26*3 is the scattering 
angle), the Born approximation is 



a = — CTt 7^— k costly 
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Aa = -2ibdZ^/^mn/me = 2 x lO'^d 
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which leads to an spin rotation, on interference with the 
Si nuclear scattering amplitude (oq = 4 x 10""'^^ cm) 



where 7 is the neutron magnetic moment (-3 Hz/mG). 
Thus, the integral is identical to the EDM case, with a 
different multipliative constant, and 



A0 = — = 4.7 X lO^^d/cm 



(10) 



Aa = —2iTrk cos 9s fJ-f{(3q) 



(18) 



implying that for d = 5 x 10 ^^ cm, a Bragg refiection 
from an Si crystal would give a rotation of 2 x 10^^" rad. 



or a spin rotation about x of 

A0 « 4 X IQ-^'^cm/A X lO^^^cm = 10"^ rad (19) 
per Bragg reflection from an Si crystal. 



V. SEMI-CLASSICAL MODEL 

If we assume there is no electron cloud around the nu- 
cleus, in Eq. (7) X = 1; ^-i^d ctedm oc 1/g, which is equiv- 
alent to the high-momentum limit. We can estimate the 
EDM effect by taking a classical trajectory with impact 
parameter b relative to the nucleus. The time-integrated- 
electric-field-induced phase shift, assuming a spin in the 
the z direction (propagation in x direction) is given by 
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which is essentially the same as before, in the high mo- 
mentum limit, if we take b = a„,,„. 




FIG. 1. Dimensionless potential integral for the Born ap- 
proximation electric field interaction. 



VI. DISCUSSION 

Comparing Eqs. (10) and (19), we see that the mo- 
tional field spin rotation is on the order of 10^ times 
larger than that due to an EDM with a magnitude that 
would be of interest in an improved experiment. Unfor- 
tunately, the effects cannot be switched on and off as in 
the case of the more conventional experiments based on 
spin precession in an applied electric field. Although the 
two scattering effects are proportional to Gx and Uz re- 
spectively, discrimination between the effects relies on an 
absolute determination of the polarization and scattering 
axes. One can also be concerned with the normal nuclear 
parity violation, which can combine with a misalignment 
to produce effects that mimic T- violation. This and other 
issues relevant for a realistic EDM scattering experiment 
are similar to those relating to a study of time reversal 
violating effects in slow neutron transmission through po- 
larized matter for which the issues have been addressed 
in some detail; in particular, the constraints on near- 
perfect field and polarization alignment, and inability 
to discriminate effects due to misalignments, have been 
emphasized. |@] The scattering angles constraints in a 
Bragg scattering EDM experiment are analogous to the 
constraints on the sample polarization axis in a neutron 
transmission experiment as discussed in [Q. Given the 
constraints (e.g., scattering angle and polarization align- 
ment to 10~^ radian absolute accuracy which requires 
10^^ neutron counts to measure experimentally) achiev- 
ing any significant increase in the limit for the neutron 
EDM would seem a daunting task. 



[2] R. Golub and J.M. Pendlebury, Contemp. Phys. 13, 519 
(1972). 

[3] I.B. Khriplovich and S.K. Lamoreaux, CP Violation with- 
out Strangeness, (Springer- Verlag, Heidelberg, 1997). Sec. 
4.3.3. 

[4] T. Dombeck, Seminar, Los Alamos National Laboratory, 
1999. 

[5] Julian Schwinger, Phys. Rev. 73, 407 (1948). 

[6] C.G. ShuU, Phys. Rev. Lett. 10, 297 (1963). 

[7] L.D. Landau and E.M. Lifshitz, Quantum Mechan- 
ics(Pergamon, Oxford, 1977). 

[8] S.K. Lamoreaux and R. Golub, Phys. Rev. D 50, 5632 
(1994). 



[1] C.G. Shull and R. Nathans, Phys. Rev. Lett. 19, 384 
(1967). 



